EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 



CERN-PH-EP-2013-039 
LHCb-PAPER-2013-003 
29 April 2013 



Search for direct CP violation in 



A search for direct CP violation in D° — > h~h + (where h = K or ir) is presented 
using data corresponding to an integrated luminosity of l.Ofb -1 collected in 2011 
by LHCb in pp collisions at a centre-of-mass energy of 7TeV. The analysis uses 
D° mesons produced in inclusive semileptonic 6-hadron decays to the D°[iX final 
state, where the charge of the accompanying muon is used to tag the flavour of the 
D° meson. The difference in the CP-violating asymmetries between the two decay 
channels is measured to be 

AA C p = A C p{K-K + ) - A C p{ir~ir + ) = (0.49 ± 0.30 (stat) ± 0.14 (syst))% . 

This result does not confirm the evidence for direct CP violation in the charm sector 
reported in other analyses. 




> — >> h modes using 
semileptonic B decays 



The LHCb collaboration!!] 



Abstract 



Published in Phys. Lett. B 



© CERN on behalf of the LHCb collaboration, license CC-BY-3.0 



t Authors are listed on the following pages. 



ii 



LHCb collaboration 



R. Aaij 40 , C. Abellan Beteta 35 -™, B. Adeva 36 , M. Adinolfi 45 , C. Adrover 6 , A. Affolder 51 , 

Z. Ajaltouni 5 , J. Albrecht 9 , F. Alessio 37 , M. Alexander 50 , S. Ali 40 , G. Alkhazov 29 , 

P. Alvarez Cartelle 36 , A.A. Alves Jr 24,37 , S. Amato 2 , S. Amerio 21 , Y. Amhis 7 , L. Anderlini 17 ^ , 

J. Anderson 39 , R. Andreassen 59 , R.B. Appleby 53 , O. Aquines Gutierrez 10 , F. Archilli 18 , 

A. Artamonov 34 , M. Artuso 56 , E. Aslanides 6 , G. Auriemma 24 '™ - , S. Bachmann 11 , J.J. Back 47 , 

C. Baesso 57 , V. Balagura 30 , W. Baldini 16 , R.J. Barlow 53 , C. Barschel 37 , S. Barsuk 7 , 

W. Barter 46 , Th. Bauer 40 , A. Bay 38 , J. Beddow 50 , F. Bedeschi 22 , I. Bediaga 1 , S. Belogurov 30 , 
K. Belous 34 , I. Belyaev 30 , E. Ben-Haim 8 , M. Benayoun 8 , G. Bencivenni 18 , S. Benson 49 , 
J. Benton 45 , A. Berezhnoy 31 , R. Bernet 39 , M.-O. Bettler 46 , M. van Beuzekom 40 , A. Bien 11 , 
S. Bifani 12 , T. Bird 53 , A. Bizzeti 17 ' fe , P.M. Bj0rnstad 53 , T. Blake 37 , F. Blanc 38 , J. Blouw 11 , 
S. Blusk 56 , V. Bocci 24 , A. Bondar 33 , N. Bondar 29 , W. Bonivento 15 , S. Borghi 53 , A. Borgia 56 , 
T.J.V. Bowcock 51 , E. Bowen 39 , C. Bozzi 16 , T. Brambach 9 , J. van den Brand 41 , J. Bressieux 38 , 

D. Brett 53 , M. Britsch 10 , T. Britton 56 , N.H. Brook 45 , H. Brown 51 , I. Burducea 28 , A. Bursche 39 , 
G. Busetto 21 ' 9 , J. Buytaert 37 , S. Cadeddu 15 , O. Callot 7 , M. Calvi 20 '^', M. Calvo Gomez 35 '™, 

A. Camboni 35 , P. Campana 18 ' 37 , D. Campora Perez 37 , A. Carbone 14 ' c , G. Carboni 23 ' fc , 

R. Cardinale 19 '*, A. Cardini 15 , H. Carranza-Mejia 49 , L. Carson 52 , K. Carvalho Akiba 2 , 

G. Casse 51 , M. Cattaneo 37 , Ch. Cauet 9 , M. Charles 54 , Ph. Charpentier 37 , P. Chen 3 ' 38 , 

N. Chiapolini 39 , M. Chrzaszcz 25 , K. Ciba 37 , X. Cid Vidal 37 , G. Ciezarek 52 , P.E.L. Clarke 49 , 

M. Clemencic 37 , H.V. Cliff 46 , J. Closier 37 , C. Coca 28 , V. Coco 40 , J. Cogan 6 , E. Cogneras 5 , 

P. Collins 37 , A. Comerma-Montells 35 , A. Contu 15 , A. Cook 45 , M. Coombes 45 , S. Coquereau 8 , 

G. Corti 37 , B. Couturier 37 , G.A. Cowan 49 , D. Craik 47 , S. Cunliffe 52 , R. Currie 49 , 

C. D'Ambrosio 37 , P. David 8 , P.N.Y. David 40 , I. De Bonis 4 , K. De Bruyn 40 , S. De Capua 53 , 

M. De Cian 39 , J.M. De Miranda 1 , L. De Paula 2 , W. De Silva 59 , P. De Simone 18 , D. Decamp 4 , 

M. Deckenhoff 9 , L. Del Buono 8 , D. Derkach 14 , O. Deschamps 5 , F. Dettori 41 , H. Dijkstra 37 , 

M. Dogaru 28 , S. Donleavy 51 , F. Dordei 11 , A. Dosil Suarez 36 , D. Dossett 47 , A. Dovbnya 42 , 

F. Dupertuis 38 , R. Dzhelyadin 34 , A. Dziurda 25 , A. Dzyuba 29 , S. Easo 48 ' 37 , U. Egede 52 , 

V. Egorychev 30 , S. Eidelman 33 , D. van Eijk 40 , S. Eisenhardt 49 , U. Eitschberger 9 , R. Ekelhof 9 , 
L. Eklund 50 ' 37 , I. El Rifai 5 , Ch. Elsasser 39 , D. Elsby 44 , A. Falabella 14 - 6 , C. Farber 11 , 

G. Fardell 49 , C. Farinelli 40 , S. Farry 12 , V. Fave 38 , D. Ferguson 49 , V. Fernandez Albor 36 , 

F. Ferreira Rodrigues 1 , M. Ferro-Luzzi 37 , S. Filippov 32 , M. Fiore 16 , C. Fitzpatrick 37 , 

M. Fontana 10 , F. Fontanelli 19 '*, R. Forty 37 , O. Francisco 2 , M. Frank 37 , C. Frei 37 , M. Frosini 17 '-/' , 
S. Furcas 20 , E. Furfaro 23 , A. Gallas Torreira 36 , D. Galli 14 ' c , M. Gandelman 2 , P. Gandini 56 , 
Y. Gao 3 , J. Garofoli 56 , P. Garosi 53 , J. Garra Tico 46 , L. Garrido 35 , C. Gaspar 37 , R. Gauld 54 , 

E. Gersabeck 11 , M. Gersabeck 53 , T. Gershon 47 ' 37 , Ph. Ghez 4 , V. Gibson 46 , V.V. Gligorov 37 , 
C. G6bel 57 , D. Golubkov 30 , A. Golutvin 52 ' 30 ' 37 , A. Gomes 2 , H. Gordon 54 , 

M. Grabalosa Gandara 5 , R. Graciani Diaz 35 , L.A. Granado Cardoso 37 , E. Grauges 35 , 

G. Graziani 17 , A. Grecu 28 , E. Greening 54 , S. Gregson 46 , O. Grunberg 58 , B. Gui 56 , 
E. Gushchin 32 , Yu. Guz 34 ' 37 , T. Gys 37 , C. Hadjivasiliou 56 , G. Haefeli 38 , C. Haen 37 , 
S.C. Haines 46 , S. Hall 52 , T. Hampson 45 , S. Hansmann-Menzemer 11 , N. Harnew , 

S.T. Harnew 45 , J. Harrison 53 , T. Hartmann 58 , J. He 37 , V. Heijne 40 , K. Hennessy 51 , P. Henrard 5 , 
J.A. Hernando Morata 36 , E. van Herwijnen 37 , E. Hicks 51 , D. Hill 54 , M. Hoballah 5 , 

C. Hombach 53 , P. Hopchev 4 , W. Hulsbergen 40 , P. Hunt 54 , T. Huse 51 , N. Hussain 54 , 

D. Hutchcroft 51 , D. Hynds 50 , V. Iakovenko 43 , M. Idzik 26 , P. Ilten 12 , R. Jacobsson 37 , A. Jaeger 11 , 

E. Jans 40 , P. Jaton 38 , F. Jing 3 , M. John 54 , D. Johnson 54 , C.R. Jones 46 , B. Jost 37 , M. Kaballo 9 , 



iii 



S. Kandybei 42 , M. Karacson 37 , T.M. Karbach 37 , I.R. Kenyon 44 , U. Kerzel 37 , T. Ketel 41 , 

A. Keune 38 , B. Khanji 20 , O. Kochebina 7 , I. Komarov 38 , R.F. Koopman 41 , P. Koppenburg 40 , 
M. Korolev 31 , A. Kozlinskiy 40 , L. Kravchuk 32 , K. Kreplin 11 , M. Kreps 47 , G. Krocker 11 , 

P. Krokovny 33 , F. Kruse 9 , M. Kucharczyk 20 ' 25 '- 5 ; , V. Kudryavtsev 33 , T. Kvaratskheliya 30 ' 37 , 
V.N. La Thi 38 , D. Lacarrere 37 , G. Lafferty 53 , A. Lai 15 , D. Lambert 49 , R.W. Lambert 41 , 
E. Lanciotti 37 , G. Lanfranchi 18,37 , C. Langenbruch 37 , T. Latham 47 , C. Lazzeroni 44 , R. Le Gac 6 , 
J. van Leerdam 40 , J. -P. Lees 4 , R. Lefevre 5 , A. Leflat 31 , J. Lefrancois 7 , S. Leo 22 , O. Leroy 6 , 

B. Leverington 11 , Y. Li 3 , L. Li Gioi 5 , M. Liles 51 , R. Lindner 37 , C. Linn 11 , B. Liu 3 , G. Liu 37 , 
S. Lohn 37 , I. Longstaff 50 , J.H. Lopes 2 , E. Lopez Asamar 35 , N. Lopez-March 38 , H. Lu 3 , 

D. Lucchesi 21,9 , J. Luisier 38 , H. Luo 49 , F. Machefert 7 , I.V. Machikhiliyan 4,30 , F. Maciuc 28 , 
O. Maev 29 ' 37 , S. Malde 54 , G. Manca 15 ' d , G. Mancinelli 6 , U. Marconi 14 , R. Marki 38 , J. Marks 11 , 
G. Martellotti 24 , A. Martens 8 , L. Martin 54 , A. Martin Sanchez 7 , M. Martinelli 40 , 
D. Martinez Santos 41 , D. Martins Tostes 2 , A. Massafferri 1 , R. Matev 37 , Z. Mathe 37 , 

C. Matteuzzi 20 , E. Maurice 6 , A. Mazurov 16 ' 32 ' 37 ' 6 , J. McCarthy 44 , R. McNulty 12 , A. Mcnab 53 , 
B. Meadows 59 ' 54 , F. Meier 9 , M. Meissner 11 , M. Merk 40 , D.A. Milanes 8 , M.-N. Minard 4 , 

J. Molina Rodriguez 57 , S. Monteil 5 , D. Moran 53 , P. Morawski 25 , M.J. Morello 22 ' 5 , 

R. Mountain 56 , I. Mous 40 , F. Muheim 49 , K. Midler 39 , R. Muresan 28 , B. Muryn 26 , B. Muster 38 , 

P. Naik 45 , T. Nakada 38 , R. Nandakumar 48 , I. Nasteva 1 , M. Needham 49 , N. Neufeld 37 , 

A.D. Nguyen 38 , T.D. Nguyen 38 , C. Nguyen-Mau 38 ' p , M. Nicol 7 , V. Niess 5 , R. Niet 9 , N. Nikitin 31 , 

T. Nikodem 11 , A. Nomerotski 54 , A. Novoselov 34 , A. Oblakowska-Mucha 26 , V. Obraztsov 34 , 

S. Oggero 40 , S. Ogilvy 50 , O. Okhrimenko 43 , R. 01deman 15 < d , M. Orlandea 28 , 

J.M. Otalora Goicochea 2 , P. Owen 52 , A. Oyanguren 35 >°, B.K. Pal 56 , A. Palano 13 - b , 

M. Palutan 18 , J. Panman 37 , A. Papanestis 48 , M. Pappagallo 50 , C. Parkes 53 , C.J. Parkinson 52 , 

G. Passaleva 17 , G.D. Patel 51 , M. Patel 52 , G.N. Patrick 48 , C. Patrignani 19 '\ C. Pavel-Nicorescu 28 , 

A. Pazos Alvarez 36 , A. Pellegrino 40 , G. Penso 24 '', M. Pepe Altarelli 37 , S. Perazzini 14 ' c , 

D. L. Perego 20j , E. Perez Trigo 36 , A. Perez-Calero Yzquierdo 35 , P. Perret 5 , M. Perrin-Terrin 6 , 

G. Pessina 20 , K. Petridis 52 , A. Petrolini 19 '\ A. Phan 56 , E. Picatoste Olloqui 35 , B. Pietrzyk 4 , 
T. Pilar 47 , D. Pinci 24 , S. Playfer 49 , M. Plo Casasus 36 , F. Polci 8 , G. Polok 25 , A. Poluektov 47 - 33 , 

E. Polycarpo 2 , D. Popov 10 , B. Popovici 28 , C. Potterat 35 , A. Powell 54 , J. Prisciandaro 38 , 
V. Pugatch 43 , A. Puig Navarro 38 , G. Punzi 22 ' r , W. Qian 4 , J.H. Rademacker 45 , 

B. Rakotomiaramanana 38 , M.S. Rangel 2 , I. Raniuk 42 , N. Rauschmayr 37 , G. Raven 41 , 

S. Redford 54 , M.M. Reid 47 , A.C. dos Reis 1 , S. Ricciardi 48 , A. Richards 52 , K. Rinnert 51 , 
V. Rives Molina 35 , D.A. Roa Romero 5 , P. Robbe 7 , E. Rodrigues 53 , P. Rodriguez Perez 36 , 
S. Roiser 37 , V. Romanovsky 34 , A. Romero Vidal 36 , J. Rouvinet 38 , T. Ruf 37 , F. Ruffini 22 , 

H. Ruiz 35 , P. Ruiz Vails 35 - , G. Sabatino 24 ' fc , J.J. Saborido Silva 36 , N. Sagidova 29 , P. Sail 50 , 

B. Saitta 15 ' rf , C. Salzmann 39 , B. Sanmartin Sedes 36 , M. Sannino 19 '*, R. Santacesaria 24 , 

C. Santamarina Rios 36 , E. Santovetti 23 ' fc , M. Sapunov 6 , A. Sarti 18 '', C. Satriano 24 '™, A. Satta 23 , 
M. Savrie 16 ' 6 , D. Savrina 30 ' 31 , P. Schaack 52 , M. Schiller 41 , H. Schindler 37 , M. Schlupp 9 , 

M. Schmelling 10 , B. Schmidt 37 , O. Schneider 38 , A. Schopper 37 , M.-H. Schune 7 , R. Schwemmer 37 , 

B. Sciascia 18 , A. Sciubba 24 , M. Seco 36 , A. Semennikov 30 , K. Senderowska 26 , I. Sepp 52 , 

N. Serra 39 , J. Serrano 6 , P. Seyfert 11 , M. Shapkin 34 , I. Shapoval 16 ' 42 , P. Shatalov 30 , 

Y. Shcheglov 29 , T. Shears 51 - 37 , L. Shekhtman 33 , O. Shevchenko 42 , V. Shevchenko 30 , A. Shires 52 , 

R. Silva Coutinho 47 , T. Skwarnicki 56 , N.A. Smith 51 , E. Smith 54 ' 48 , M. Smith 53 , M.D. Sokoloff 59 , 

F. J.P. Soler 50 , F. Soomro 18 , D. Souza 45 , B. Souza De Paula 2 , B. Spaan 9 , A. Sparkes 49 , 

P. Spradlin 50 , F. Stagni 37 , S. Stahl 11 , O. Steinkamp 39 , S. Stoica 28 , S. Stone 56 , B. Storaci 39 , 
M. Straticiuc 28 , U. Straumann 39 , V.K. Subbiah 37 , S. Swientek 9 , V. Syropoulos 41 , 



iv 



M. Szczekowski 27 , P. Szczypka 38 ' 37 , T. Szumlak 26 , S. T'Jampens 4 , M. Teklishyn 7 , 

E. Teodorescu 28 , F. Teubert 37 , C. Thomas 54 , E. Thomas 37 , J. van Tilburg 11 , V. Tisserand 4 , 

M. Tobin 38 , S. Tolk 41 , S. Topp-Joergensen 54 , N. Torr 54 , E. Tournefier 4,52 , S. Tourneur 38 , 

M.T. Tran 38 , M. Tresch 39 , A. Tsaregorodtsev 6 , P. Tsopelas 40 , N. Tuning 40 , M. Ubeda Garcia 37 , 

A. Ukleja 27 , D. Urner 53 , U. Uwcr 11 , V. Vagnoni 14 , G. Valenti 14 , R. Vazquez Gomez 35 , 

P. Vazquez Regueiro 36 , S. Vecchi 16 , J.J. Velthuis 45 , M. Veltri 17 ^, G. Veneziano 38 , 

M. Vesterinen 37 , B. Viaud 7 , D. Vieira 2 , X. Vilasis-Cardona 35 '™, A. Vollhardt 39 , D. Volyanskyy 10 , 

D. Voong 45 , A. Vorobyev 29 , V. Vorobyev 33 , C. VoB 58 , H. Voss 10 , R. Waldi 58 , R. Wallace 12 , 

S. Wandernoth 11 , J. Wang 56 , D.R. Ward 46 , N.K. Watson 44 , A.D. Webber 53 , D. Websdale 52 , 

M. Whitehead 47 , J. Wicht 37 , J. Wiechczynski 25 , D. Wiedner 11 , L. Wiggers 40 , G. Wilkinson 54 , 

M.P. Williams 47 ' 48 , M. Williams 55 , F.F. Wilson 48 , J. Wishahi 9 , M. Witek 25 , S.A. Wotton 46 , 

S. Wright 46 , S. Wu 3 , K. Wyllie 37 , Y. Xie 49 ' 37 , F. Xing 54 , Z. Xing 56 , Z. Yang 3 , R. Young 49 , 

X. Yuan 3 , O. Yushchenko 34 , M. Zangoli 14 , M. Zavertyaev 10 ' a , F. Zhang 3 , L. Zhang 56 , 

W.C. Zhang 12 , Y. Zhang 3 , A. Zhelezov 11 , A. Zhokhov 30 , L. Zhong 3 , A. Zvyagin 37 . 

1 Centro Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 

2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 

3 Center for High Energy Physics, Tsinghua University, Beijing, China 
4 LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

5 Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
6 CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
7 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

8 LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 

9 Fakultat Physik, Technische Universitat Dortmund, Dortmund, Germany 

10 Max-Planck-Institut fur Kernphysik (MPIK), Heidelberg, Germany 

11 Physikalisches Institut, Ruprecht-Karls-Universitdt Heidelberg, Heidelberg, Germany 

12 School of Physics, University College Dublin, Dublin, Ireland 
13 Sezione INFN di Bari, Bari, Italy 

14 Sezione INFN di Bologna, Bologna, Italy 
15 Sezione INFN di Cagliari, Cagliari, Italy 
16 Sezione INFN di Ferrara, Ferrara, Italy 

17 Sezione INFN di Firenze, Firenze, Italy 

18 Laboratori Nazionali dell 'INFN di Frascati, Frascati, Italy 

19 Sezione INFN di Genova, Genova, Italy 

20 Sezione INFN di Milano Bicocca, Milano, Italy 

21 Sezione INFN di Padova, Padova, Italy 

22 Sezione INFN di Pisa, Pisa, Italy 

23 Sezione INFN di Roma Tor Vergata, Roma, Italy 
24 Sezione INFN di Roma La Sapienza, Roma, Italy 

25 Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland 
26 AGH University of Science and Technology, Krakow, Poland 
27 National Center for Nuclear Research (NCBJ), Warsaw, Poland 

2S Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

29 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

30 Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

31 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 

32 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 

33 Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 

34 Institute for High Energy Physics (IHEP), Protvino, Russia 

35 Universitat de Barcelona, Barcelona, Spain 



V 



36 Universidad de Santiago de Compostela, Santiago de Compostela, Spain 
37 European Organization for Nuclear Research ( CERN), Geneva, Switzerland 
3S Ecole Polytechnique Federate de Lausanne (EPFL), Lausanne, Switzerland 
39 Physik-Institut, Universitat Zurich, Zurich, Switzerland 

40 Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 

41 Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The 

Netherlands 

42 NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 
^Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 
44 University of Birmingham, Birmingham, United Kingdom 
45 H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

46 Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 

47 Department of Physics, University of Warwick, Coventry, United Kingdom 
48 STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

49 School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
50 School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 

51 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 

52 Imperial College London, London, United Kingdom 

53 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 
54 Department of Physics, University of Oxford, Oxford, United Kingdom 
55 Massachusetts Institute of Technology, Cambridge, MA, United States 
56 Syracuse University, Syracuse, NY, United States 

57 Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 
5S Institut fur Physik, Universitat Rostock, Rostock, Germany, associated to 11 
59 University of Cincinnati, Cincinnati, OH, United States, associated to 56 

a P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 

h Universitd di Bari, Bari, Italy 

c Universitd di Bologna, Bologna, Italy 

d Universitd di Cagliari, Cagliari, Italy 

e Universitd di Ferrara, Ferrara, Italy 

f Universitd di Firenze, Firenze, Italy 

9 Universitd di Urbino, Urbino, Italy 

Universitd di Modena e Reggio Emilia, Modena, Italy 
1 Universitd di Genova, Genova, Italy 
3 Universitd di Milano Bicocca, Milano, Italy 

Universitd di Roma Tor Vergata, Roma, Italy 
1 Universitd di Roma La Sapienza, Roma, Italy 
m Universitd delta Basilicata, Potenza, Italy 

n LIFAELS, La Salle, Universitat Ramon Hull, Barcelona, Spain 

°IFIC, Universitat de Valencia- C SIC, Valencia, Spain 

p Hanoi University of Science, Hanoi, Viet Nam 

q Universitd di Padova, Padova, Italy 

r Universitd di Pisa, Pisa, Italy 

s Scuola Normale Superiore, Pisa, Italy 



vi 



1 Introduction 



The combined symmetry of charge conjugation and parity (CP) is broken in the weak 
interaction of the Standard Model by a single phase in the Cabibbo-Kobayashi-Maskawa 
matrix |1] |2j|. Physics beyond the Standard Model may reveal itself in the form of 
additional sources of CP violation. In both the K° and B° systems CP violation has been 
unambiguously observed, and is in agreement with the Standard Model predictions. In 
contrast, CP violation in the charm sector has yet to be established. The amount of CP 
violation in charm decays was generally expected to be much smaller than the 1% level 
in the Standard Model [21 E]. The LHCb collaboration, however, reported evidence with 
3.5 standard deviations significance for direct CP violation in two-body, singly-Cabibbo- 
suppressed D° decays [5]. The difference in CP asymmetries between D°— > K~K + and 
D° — > 7T-7T+ decays was found to be AA CP = (-0.82 ± 0.21 (stat) ± 0.11 (syst))%. This 
result sparked a theoretical debate on whether or not this could be accommodated within 
the Standard Model. For a comprehensive review see Ref. [6]. 

After the LHCb paper, the CDF and Belle collaborations presented measurements 
of AA CP = (-0.62 ± 0.21 (stat) ± 0.10 (syst))% and AA CP = (-0.87 ± 0.41 (stat) ± 
0.06 (syst))% [8], respectively. These numbers are included in the average from the Heavy 
Flavor Averaging Group (HFAG) [9], together with a previous measurement pU] from the 
BaBar collaboration, yielding a world average of the difference in direct CP violation of 
Aagj, = (-0.68 ±0.15)% J] 

In all previous results D* + — > D°tt + decays^] have been used as the source of the D° 
sample, and the emitted pion was used to determine the flavour of the neutral D meson 
(i.e., whether it is D° or D°). In this paper a measurement of AAcp is presented using D° 
mesons produced in semileptonic 6-hadron decays where the flavour of the neutral D meson 
is tagged by the accompanying charged lepton. This approach provides an independent 
determination of AAqp- 



2 Method and formalism 

The measured (raw) asymmetry for a D° decay to a CP eigenstate / is defined as 

_ N(D° -> f) - N(D° -> f) 
raw N(D° -> /) + N(D° -> /) ' 1 ' 

where N denotes the observed yield for the given decay. The initial flavour of the neutral 
D meson is tagged by the charge of the accompanying muon in the semileptonic 6-hadron 
(B) decay to the D/j,X final state. A positive muon is associated with a D° meson, and a 
negative muon with a D° meson. The X denotes any other particle(s) produced in the 
semileptonic B decay, which are not reconstructed (e.g., the neutrino). 



1 The relation between AAcp and Aa^lp is explained in Sect. 6j 
2 The inclusion of charge-conjugated modes is implied throug 



rout this paper, unless explicitly stated 



otherwise. 
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The raw asymmetry can be written in terms of the D° decay rate, T, the muon detection 
efficiency, e, and the D° production rate in semileptonic 6-hadron decays, V, as 

= T(D°)e(j J .-)V(D ) - T(D°)e(p+)V(D°) 
mw T(D°)e(fj,-)V(D°) + T(D°)e(fj,+)V(D°) ' U 

Defining the CP asymmetry as A CP = (T(D°) - T(D°))/(T(D ) + T(D )), the muon 
detection asymmetry as A^ D = (e(fi")—e(fi + j)/(e(fi + )+e(fi~)), and the effective production 
asymmetry as Ap = (V(D°) — V (D )) / (V (D°) + V(D )), the raw asymmetry can be 
written to first order as 

A vaw « A CP + Ap + Ap . (3) 

The effective production asymmetry is due to different production rates of b- and &-hadrons 
and also includes any effect due to semileptonic asymmetries in neutral B mesons. As the 
detection and production asymmetries are of order 1%, the approximation in Eq. ([3]) is 
valid up to corrections of order 10~ 6 . Both detection and production asymmetries differ 
from those in the analyses using D*^ decays, where the D*^ 1 mesons are produced directly 
in the primary pp interaction. In these "prompt" decays a possible detection asymmetry 
enters through the reconstruction of the tagging pion, and the production asymmetry is 
that of the prompt D*^ 1 mesons. 

By taking the difference between the raw asymmetries measured in the K~K + 
and D° — > 7r~7r + decays the detection and production asymmetries cancel, giving a robust 
measurement of the CP asymmetry difference 

AA CP = A iaw (K-K + ) - A raw (7r-7r+) « A CP (K- K + ) - A CP {ix-^) . (4) 

Since the detection and the production depend on the kinematics of the process under 
study, the cancellation is only complete when the kinematic distributions of the muon and 
6-hadron are the same for both D° — > K~K + and D° — > 7r~7r + . A weighting procedure is 
used to improve the cancellation by equalising the kinematic distributions. 



3 Detector and simulation 

The LHCb detector [If J is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The polarity of the magnet is reversed 
repeatedly during data taking, which causes all detection asymmetries that are induced by 
the left-right separation of charged particles to change sign. The combined tracking system 
has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and impact parameter resolution of 20 |xm for tracks with high transverse momentum. 
Charged hadrons are identified using two ring- imaging Cherenkov detectors [12] . Muons 
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are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. The trigger jT3] consists of a hardware stage, based on information from the 
calorimeter and muon systems, followed by a software stage, which applies a full event 
reconstruction. 

In the simulation, pp collisions are generated using Pythia 6.4 [T4] with a specific 
LHCb configuration [15]. Decays of hadronic particles are described by EvtGen [16] 
in which final state radiation is generated using Photos [T7]. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 
toolkit [18] as described in Ref. [15] . The B + and B° mesons in the simulated events are 
forced to decay semileptonically using a cocktail of decay modes, including those that 
involve excited D states and r leptons, that lead to final states with a D° meson and a 
muon. 

4 Data set and selection 

This analysis uses the LHCb 2011 data set, corresponding to an integrated luminosity of 
1.0 fb~ , of which 0.4fb _1 is taken with the magnet field pointing up and 0.6 fb -1 with the 
magnet field pointing down. The measurement of AAcp is performed separately for the 
two field polarities. The final value for AAcp is obtained by taking the arithmetic mean of 
the two results to reduce as much as possible any residual effect of the detection asymmetry. 
To minimise potential trigger biases the candidates are required to be accepted by specific 
trigger decisions. About 87% of the candidates in the final selection are triggered at the 
hardware stage by the muon system only, about 3% by the hadronic calorimeter only and 
about 10% by both. The muon trigger requires the muon transverse momentum, p^, to 
be greater than 1.48GeV/c. The effect of a charge-dependent shift in the pt estimate in 
this trigger is corrected, which requires tightening the muon transverse momentum cut, 
as measured by the hardware trigger, to > 1.64GeV/c. In the software trigger the 
candidates are selected by either a single muon trigger or by a topological trigger, which 
selects combinations of a muon with one or two additional tracks that are consistent with 
the topological signature of 6-hadron decays. At this level, 5% of the candidates in the 
final selection are selected by the single muon trigger only, 79% by the topological trigger 
only, and 16% by both. 

In order to suppress backgrounds, the y 2 P er degree of freedom of the track fit is 
required to be smaller than 4 for the kaons and pions and smaller than 5 for the muon. 
Furthermore, the y 2 per degree of freedom of each of the 6-hadron and D° decay vertex fits 
is required to be smaller than 6, and the impact parameter y 2 (defined as the difference 
between the y 2 of the primary vertex formed with and without the considered tracks) is 
required to be larger than 9 for all three tracks. The significance of the distance between 
the primary vertex and the D° decay vertex is required to be above 10. The momentum 
and transverse momentum of the muon are required to be above 3 GeV/c and 1.2 GeV/cJ^] 
and the momentum and transverse momentum of the D° daughters above 2 GeV/c and 

3 This cut affects mainly the candidates triggered by the hadronic calorimeter at the hardware level. 
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0.3GeV/c. The D° transverse momentum must be above 0.5GeV/c and the scalar pr sum 
of its daughters above 1.4GeV/c. The invariant mass of the D°-muon combination is 
required to be between 2.5 and 5.0GeV/c 2 to suppress background. The upper bound 
removes three-body final state 6-hadron decays. The reconstructed decay time of the 
D° meson (measured from the fe-hadron decay vertex) is required to be positive. The 
requirement on the muon impact parameter reduces the contribution from D° mesons 
produced directly in the pp collision to below 3%. Requirements on the D° decay topology 
are minimal in order to keep the lifetime acceptance similar for the D° — > K~ K + and 
D° — > 7r~7r + modes. 

A potentially significant background from B — > J/ip X decays is suppressed by removing 
candidates where the invariant mass of the muon and the oppositely-charged D° daughter 
is within three times the mass resolution from the J/ip or ip(2S) mass and the D° 
daughter passes muon identification requirements. Reflections from Cabibbo-favoured 
K~tt + decays are observed in the mass regions below and above the signal peaks in 
the D° — > 7r~7r + and D° — > K~K + samples, respectively. Information from the relevant 
detectors in LHCb is combined into differences between the logarithms of the particle 
identification likelihoods under different mass hypotheses (DLL). The selected kaons are 
required to have DLL^ = \nC K — ln£ w > 10 and the selected pions are required to have 
DLL^^ < —2. The D° — > K~n + mode is used as a control channel and is selected with 
the same requirements as the two decay modes of interest. 

5 Determination of the asymmetries 

The invariant mass distributions for the muon-tagged D° candidates are shown in Fig. [I] To 
determine the numbers of signal candidates after selection, a binned maximum likelihood 
fit to each of these distributions is performed. The signal is modelled by the sum of 
two Gaussian functions with common means, but different widths. The combinatorial 
background is described by an exponential shape. For the tt~7t + invariant mass distribution 
the fit is performed in the range between 1795 — 1940MeV/c 2 and a Gaussian distribution 
is used to model the tail of the reflection from D° — > K~ix + decays. For the K~K + 
invariant mass distribution the fit range is restricted to 1810 — 1920MeV/c 2 such that 
the contamination from the D° — > K~n + reflection and from partially reconstructed 
D° — > K~K + 7T° and D + — > K~K + it + decays is negligible. The total number of signal 
candidates determined from the fit is (558.9 ± 0.9) x 10 3 for D° — > K~K + decays and 
(221.6 ± 0.8) x 10 3 for D°— > 7^+ decays. 

The raw asymmetries are determined with simultaneous binned likelihood fits to the 
D° mass distributions for positive and negative muon tags where the shape parameters for 
the signal and the D -^ K~ir + reflection are required to be the same. The background 
shape can vary independently for positive and negative muon tags. Table [l] lists the 
raw asymmetries for both modes, and for the D° — > K~tt + control mode. An additional 
asymmetry in the D°— > K~tt + mode originates from the different cross-sections in matter 
for positive and negative kaons. It can be seen that the asymmetry in this mode is consistent 
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Figure 1: Invariant mass distributions for (a, c) I? —)- K~ K + and (b, d) D — > tt~tt + muon- 
tagged candidates for the two magnet polarities. The result of the fit is overlaid, showing the 
contribution from signal, combinatorial background and D° — > K~n + reflection. Underneath 
each plot the pull in each mass bin is shown. 

for the two magnetic field polarities, which indicates that the detection asymmetry related 
to the magnetic field is at most C(10~ 3 ). 

Table 1: Unweighted raw asymmetries (in %) for the D°— > ir~7r + , D° — > K~ K + and D° — > K~tt + 
decays for the two magnet polarities. The mean value is the arithmetic average over the two 
polarities. The uncertainties are statistical only. 





Magnet up 


Magnet down 


Mean 


unweighted ^ j£ — ^ 

A unweighted / - +\ 
raw V / 


-0.33±0.23 
-1.18±0.40 


-0.22±0.19 
-0.35±0.34 


-0.28±0.15 
-0.77±0.26 




0.85±0.46 


0.13±0.39 


0.49±0.30 


^unweighted 


-1.64±0.10 


-1.60±0.08 


-1.62±0.06 
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5.1 Differences in kinematic distributions 



Since the detection and production asymmetries may have kinematic dependences, the 
cancellation in Eq. (|4j) is only valid if the kinematic distributions of the muon and b- 
hadron are similar for both D°—> K~K + and D° — > tt~7t + decays. After the trigger and 
selection requirements the kinematic distributions for the two decay modes are, however, 
slightly different. Although the selection is largely the same, the particle identification 
requirements introduce differences in the momentum distributions. In addition, due to 
the difference in available phase space, the pions in D° — > 7t~7t + decays have a harder 
momentum spectrum compared to the kaons in D°— >■ K~K + decays. The muon trigger 
and selection requirements are identical. Nevertheless, the D° meson and the muon are 
kinematically correlated since they originate from the same decay, causing also the muon 
kinematic distributions to be different for the two decay modes. Figure [2] shows the px and 
pseudorapidity 77 distributions for the D° meson and the muon. The background has been 
statistically subtracted using the sPlot method [20j. In order to obtain the same kinematic 
distributions for both decays, the D° candidates are given a weight depending on their pr 
and t] values. The weights are obtained from a comparison of the background-subtracted 
distributions and are applied to either D°— > K~K + or D° — > 7r~7r + candidates depending 
on which has most events in the given kinematic bin. Figure [3] shows the weighted 
kinematic distributions for both decay modes. Whereas the weights are determined purely 
on the basis of the D° and rj distributions, after the weighting, the muon distributions 
are also in excellent agreement. The raw asymmetries after the weighting procedure for 
the D°— >■ K~K + and _D°— > tt--k + modes are given in Table [2] There are minor changes 
in the values of the raw asymmetries and AAqp with respect to the unweighted results, 
showing that the effect of the difference in kinematic distributions is small. 

Table 2: Weighted raw asymmetries (in %) for the D°— > tt~it + and K~ K + decays for the 

two magnet polarities. The mean value is the arithmetic average over the two polarities. The 
uncertainties are statistical only. 





Magnet up 


Magnet down 


Mean 


A raw (K-K+) 


-0.39±0.23 


-0.20±0.19 


-0.29±0.15 


A-aw(7T~7r + ) 


-1.25±0.40 


-0.29±0.34 


-0.77±0.26 


AA C P 


0.86±0.46 


0.09±0.39 


0.48±0.30 



5.2 Wrong flavour tags 

In some cases the D° flavour is not tagged correctly by the muon charge due to misrecon- 
struction (e.g., a prompt D° decay can be combined with a random muon). The probability 
to tag a D° meson with a positive muon is denoted by u + and the probability to tag a D° 
meson with a negative muon by u~ . The average mistag probability is u = (u> + + ui~)/2 
and the mistag difference is Aco = co + — u~ . The raw asymmetry in Eq. d3j) is then 
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Figure 2: Kinematic distributions of the (a, c) D° meson and (b, d) muon for D°— > it~it + (black 
circles) and D° — >■ K~ K + (red squares) candidates normalised to unit area. The histograms 
show the distributions of signal candidates, after background subtraction. Underneath each plot 
the ratio of the two distributions is shown. 



modified to 

A iaw n(l-2u)(A CP + A> x D + A*)-Auj , (5) 

which makes clear that the average mistag probability dilutes the observed asymmetry, 
while any difference in the mistag probability for D° and D° gives rise to a systematic 
shift in A rsm . Assuming that the values of u and Ao; are the same for D° — > K~K + and 
D° — > 7r~7r + , the value of AAqp is then corrected as 

AA CP = (1 - 2u)- 1 (A raw (K-K + ) - A raw (7r-7r + )) . (6) 

The mistag probability is estimated from the D°— > K~tt + sample. As the D° — > K~n + 
decay is almost self-tagging the mistag probability is determined using the charge of the 
final state (either K + n~ or K~n + ). The wrongly tagged decays include a fraction of 
doubly-Cabibbo-suppressed D° — » K + tt~ and mixed D° — > D° — > K + n~ decays. This 
fraction is calculated to be (0.393 ± 0.007)% using input from Ref. [21]. After correcting 
for this fraction the average mistag probability, uJ, is found to be (0.982 ± 0.012)%, 
which means that the effect from wrong tags constitutes only a small correction on the 
observed asymmetries. This number also provides an upper bound of about 2% from any 
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Figure 3: Kinematic distributions of the (a, c) D meson and (b, d) muon for D°— >■ 7r~7r + (black 
circles) and D°— )■ K~ K + (red squares) candidates normalised to unit area after the weighting 
procedure. The histograms show the distributions of signal candidates, after background 
subtraction. Underneath each plot the ratio of the two distributions is shown. 



background from real D° decays with a random muon, which includes promptly produced 
D° decays. The difference in mistag probabilities for D° and D° mesons is found to be 
Acu = (0.006 ± 0.021)% and is neglected. 

As a cross-check the mistag probabilities are also determined from a doubly-tagged 
sample by reconstructing B — > D* + fi~X decays where the D* + decays to D°tt + and 
comparing the charge of the pion with that of the muon. The fraction of wrongly tagged 
decays is estimated from a simultaneous fit, similar to that in Ref. [22], to the distribution 
of AM = M(h~ h + -K + ) — M(h~h + ) for the full sample and for the wrongly tagged decays. 
The mistag probability in the D°— >■ K~n + sample is (0.880 ± 0.043)%, while the average 
mistag probability in the D° — > K~K + and D° — > 7r~7r + samples equals (1.00 ± 0.09)%. 
The largest difference with the result obtained from the full D° — > K~tc + sample (i.e., 
0.102%) is assigned as a systematic uncertainty in the mistag probability. The difference 
in mistag probabilities, Acu, in this cross-check is also consistent with zero. 

After the weighting and correcting for the mistag probability of (0.982 ± 0.012 (stat) ± 
0.102 (syst))%, the difference of the raw asymmetries between the two modes is found to 
be 

AA CP = (0.49 ± 0.30)% , 
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where the uncertainty is statistical only. The corresponding systematic uncertainties are 
discussed in Sect. 



6 Measurement of the average decay times 

The time-integrated asymmetry for a decay to a CP eigenstate / is defined as 

_ r(D° -> f) - t(d° -> f) 
CP t(d° ->•/) + r(£>° ->►/)' u 

where V is the decay rate for the given channel. As the reconstruction and selection 
requirements for the two decay modes are not identical, the decay time acceptance can 
be different. This introduces a difference in the contribution from direct and indirect CP 
violation for the two modes. When assuming the CP violating phase in D° oscillations, 
cf), to be universal [3], the difference between the asymmetries for D° — > K~K + and 
D° — > 7r~7r + can be written in terms of direct and indirect CP violation as |23] 




Aa^ r P I 1 + y^- cos 1 + ( a££ + a% P y cos <p ) — ^ . (8) 



A(t) 



In this equation the indirect CP violation is a^p = — (A m /2)y cos<f) + xsin0, x and y are 
the D° mixing parameters, A m represents the CP violation from mixing, r is the average 
D° lifetime, Aa^j and a^ p are the direct CP violation difference and average of the two 
decay modes, and A(t) and (t) are the difference and average of the two mean decay 
times. Under SU(3) flavour symmetry, the direct asymmetries in the individual modes 
are expected to have opposite sign and therefore add constructively in the difference. 
Furthermore, since y is of order 1%, (t)/r is 0(1) and A(t)/r is close to zero, AAcp is 
essentially equal to the difference in direct CP violation, Aa^p. While y and cos0 can be 
obtained from the HFAG averages [9] , in order to interpret A Aqp in terms of direct and 
indirect CP violation, the mean decay time (t) in each channel needs to be measured. 

The determination of the mean decay time is performed through a fit to the decay time 
distribution of the signal candidates. Candidates with negative measured decay times are 
included in the fit to have a better handle on the acceptance and the resolution function. 
The measured decay time distribution is modelled by a decreasing exponential function, 
with mean lifetime r, convolved with a double Gaussian resolution function and multiplied 
with an acceptance function of the form 

A(t) = 1 - ae- {t/{bT))2 , (9) 

where a and b are acceptance parameters. The fit model is motivated by simulation 
studies. The values for the fraction and width of the second Gaussian and the acceptance 
parameter b are taken from the simulation and fixed in the fit. The role of the acceptance 
parametrisation is to allow a fit to the distribution such that the resolution effect can 
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Figure 4: Decay time distribution for signal candidates (solid points) with the result from the 
fit overlaid for (a) D°— > K~ K + and (b) D°— > -k~it + decays. The distribution for background 
candidates scaled to a ±34 MeV/c 2 window around the nominal D° mass is shown in the shaded 
(green) region. The distributions for signal and background candidates are obtained using the 
sPlot method. 



be removed and the true decay time, which appears in Eq. (J8J), can be evaluated. The 
observed decay time distributions with the fit result superimposed are shown in Fig. |4| 

The decay time resolutions obtained from the lifetime fit (taken as the width of the first 
Gaussian function) are 63.3±0.3fs for K~K + and 58.3±0.4fs for D -^ 7r~7r + , which 

are about 10% larger than expected from simulations. The main systematic uncertainties 
come from the uncertainty in the acceptance function and from backgrounds. Using the 
world average of the D° lifetime, t(D°) = 410.1 ± 1.5 fs, the difference and average of the 
mean decay times relative to t(D°) are found to be 

A(t)/r{D ) = 0.018 ± 0.002 (stat) ± 0.007 (syst) 
(t)/r(Z?°) = 1.062 ± 0.001 (stat) ± 0.003 (syst) , 

where the uncertainty in t(D°) is included as a systematic uncertainty. Note that (t) is 
not a measurement of the D° effective lifetime (i.e., the lifetime measured with a single 
exponential fit), since this number contains effects from the LHCb acceptance. The small 
value of A(t) implies that the measured value of AAcp is equal to the difference in direct 
CP violation, i.e., AAcp = Aa^p with negligible corrections. 

7 Systematic uncertainties 

The contributions to the systematic uncertainty on AAcp are described below. 

• Difference in b-hadron mixture. Due to the momentum requirements in the 
trigger and selection, the relative contribution from B° and B + decays (the con- 
tribution from 6-baryon and B° s decays can be neglected) can be different between 
the D°— >■ K~K + and D° — > 7r~7r + modes. In combination with a different effective 



10 



production asymmetry for candidates from B° and B + mesons (the production asym- 
metry from B° mesons is diluted due to B° mixing) this could lead to a non- vanishing 
bias in AAcp- Assuming isospin symmetry, the production cross-sections for B° and 
B + mesons are expected to be equal. Therefore, the ratio between B° and B + decays 
is primarily determined by their branching fractions to the _D°/iX final state. Using 
the inclusive branching fractions [24|, B +, ° — » D°X, the B° fraction is expected to 
be f(B°) = (37.5 ± 2.9)%. From the simulation the difference in the B° fraction due 
to the difference in selection efficiencies is found to be at maximum 1%. Further 
assuming a B + production asymmetry of 1.0% [25] and assuming no B° production 
asymmetry, the difference in the effective production asymmetry between the two 
modes is ~ 0.02%. 

Difference in B decay time acceptance. A difference between the D°— > K~ K + 

and D° — > 7t~tt + modes in the B decay time acceptance, in combination with B° 
mixing, changes the effective B production asymmetry. Its effect is estimated from 
integrating the expected B decay time distributions at different starting values, such 
that the mean lifetime ratio corresponds to the observed B decay length difference 
(~ 5%) in the two modes. Using the estimated B° fraction and assuming a 1.0% 
production asymmetry, the effect on AAcp is found to be 0.02%. 

Effect of the weighting procedure. After weighting the D° distributions in pr 
and 7], only small differences remain in the muon kinematic distributions. In order to 
estimate the systematic uncertainty from the B production and detection asymmetry 
due to residual differences in the muon kinematic distributions, an additional weight 
is applied according to the muon (pt,v) an d the azimuthal angle 0. The value of 
AA C p changes by 0.05%. 

Difference in mistag asymmetry. The difference in the mistag rate between 
positive and negative tags contributes to the measured raw asymmetry. The mistag 
difference using D° — > K~n + decays is measured to be Auj = (0.006 ± 0.021)% (see 



Sect. 5.2). In case Aw is different for D° — > K~K + and D° — > tt~7t + there can be 
a small effect from the mistag asymmetry. A systematic uncertainty of 0.02% is 
assigned, coming from the uncertainty on Acu. 

Effect of different fit models. A possible asymmetry in the background from false 
D° combinations is accounted for in the fit to the D° mass distribution. Different 
models can change the fraction between signal and background and therefore change 
the observed asymmetry. The baseline model is modified by either using a single 
Gaussian function for the signal, a single Gaussian plus a Crystal Ball function for 
the signal, a first- or second-order polynomial for the background, by leaving the 
asymmetry in the reflection free, or by modifying the fit range for D° — > tt~tt + to 
exclude the reflection peak. The largest variation changes the value of AAcp by 
0.035%. As another check, the asymmetry is determined without any fit by counting 
the number of positively- and negatively-tagged events in the signal window and 
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subtracting the corresponding numbers in the sideband windows. The sideband 
windows are defined as [/z sig — 48 MeV/ c 2 , /i sig — 34 MeV/ c 2 ] and [/i sig + 34 MeV/ c 2 , /i sig + 
48MeV/c 2 ], and the signal window as [/i sig - 14 MeV/ c 2 , /i sig + 14 MeV/ c 2 ], where /i sig 
is the mean of the signal distribution. This method changes the value of AAqp by 
0.05%, which is taken systematic uncertainty. 

• Low-lifetime background in D° — > 7r~7r + . As can be seen in Fig. |4j there 
is more background around t = in the D° — > 7r~7r + channel compared to the 
D°— > K~K + channel. If this background exhibits a non-flat or peaking structure 
this could bias the measurement of AAqp- When including the negative lifetime 
events the value of AAcp changes by 0.11%. This shift is taken as a systematic 
uncertainty. 

• background in D° — y K~K + . A non- negligible fraction of the background 
in the D° — > K~K + mode originates from partial reconstruction of A£ — > pK~7i + 
decays, where the proton is misidentified as a kaon. Most of these A£ decays are 
expected to come from semileptonic A® decays. From exclusively reconstructed A£ 
decays the shape of the background is observed to be linear in the K~K + invariant 
mass distribution. The influence of such a linear background on the fit model is 
tested by generating many pseudo-experiments. With an asymmetry in the A£ 
background of 3%, which is a conservative upper bound for the asymmetry observed 
in the exclusively reconstructed A£ decays, a small bias of 0.03% is seen in the 
measured asymmetry. This bias is taken as a systematic uncertainty. 

The systematic uncertainties are summarised in Table |3j The effects from higher-order 
corrections to Eq. ^ and of the uncertainty in the average mistag rate are found to be 
negligible. The overall systematic uncertainty on AAcp, obtained by adding the individual 
contributions in quadrature, is 0.14%. 

8 Cross-checks 

Many cross-checks have been performed to verify the stability of the result. In particular, 
the raw asymmetries and AAcp are found to be stable when applying fiducial cuts in 
the two-dimensional space of the muon momentum and its horizontal component, when 
comparing different trigger decisions and when applying tighter particle identification 
requirements on the D° daughters or on the muons. The stability of the raw asymmetries 
and AAcp is also investigated as a function of all possible reconstructed quantities, for 
instance the D° decay time, the 6-hadron flight distance, the reconstructed _D°-muon 
mass, the angle between the muon and D° daughters, and the (transverse) momenta and 
pseudorapidity of the muon and D° meson. No significant dependence is observed in 
any of these variables. For example, Fig. [5] shows AAcp and the raw asymmetries in the 
D° — > K~K + and D° — > 7t~tt + modes as a function of px and i] of the D° meson, which 
are the variables that are used in the weighting procedure. To check for a possible time 
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Table 3: Contributions to the systematic uncertainty of AAcp- 





Absolute 


Source of uncertainty 


uncertainty 


Production asymmetry: 




Difference in 6-hadron mixture 


0.02% 


Difference in B decay time acceptance 


0.02% 


Production and detection asymmetry: 




Different weighting 


0.05% 


Background from real D° mesons: 




Mistag asymmetry 


0.02% 


Background from fake D° mesons: 




D° mass fit model 


0.05% 


Low-lifetime background in D° — > 7t~tt + 


0.11% 


A£ background in D° — > K~K + 


0.03% 


Quadratic sum 


0.14% 
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dependence of the detection asymmetry the data taking period is divided into six parts of 
roughly equal integrated luminosity. The six parts are separated by periods without beam 
and changes in the magnet polarity. No significant variation of the raw asymmetries is 
observed. 



9 Conclusion 

The difference in CP asymmetries between the D° — > K~K + and D° — > tt~tt + modes is 
measured using D° mesons produced in semileptonic B decays and is found to be 

AA CP = (0.49 ± 0.30 (stat) ± 0.14 (syst))% . 
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This result takes into account the muon mistag probability and differences in the kinematic 
distributions of D° — > K~K + and D° — > 7t~tt + decays. When neglecting indirect CP 
violation the difference between this result and the previous published LHCb result using 
prompt D° decays [5] is 3.2 standard deviations, assuming that the uncertainties have a 
Gaussian distribution. The discrepancy, however, is reduced to 2.2 standard deviations 
comparing to a preliminary update of the previous result [26] . This result does not confirm 
the evidence for direct CP violation in the charm sector. 
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